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Abstract. In light of the recent finding of the narrow clustering of the geometrically-corrected 
gamma-ray energies emitted by Gamma Ray Bursts (GRBs), we investigate the possibility to use 
these sources as standard candles to probe cosmological parameters such as the matter density Q.„, 
and the cosmological constant energy density Q.^. By simulating different samples of gamma-ray 
bursts, based on recent observational results, we find that £l m (with the prior £2 m +£1a = 1) can 
be determined with accuracy ^7% with data from 300 GRBs, provided a local calibration of the 
standard candles be achieved. 



INTRODUCTION 

Recent studies have pointed out that Gamma-Ray Bursts (GRBs) may be considered as 
standard cosmological candles. The prompt y-ray energies of GRBs, after correction for 
the conical geometry of the jet, result clustered around a mean value of a few 10 50 erg 

n. 

Since the discovery that GRBs lie at cosmological distances, about 30 redshifts have 
been measured. Apart from the controversial case of GRB 980425, possibly associated 
with the nearby supernova SN1998bw (at z = 0.0085), all other redshifts are spread 
within the wide 0.17—4.5 range. Therefore GRBs could be good candles to probe 
cosmological parameters 0] 0]. 

GRBs are thought to be associated with the death of massive (and short lived) stellar 
progenitors. Therefore the rate of GRB events per unit cosmological volume should be 
a tracer of the global history of star formation. 

Hence, we have now all the information necessary to perform simulations of GRB dis- 
tributions in a given cosmological model. Universal parameters such as the matter den- 
sity fraction Q. m and the cosmological constant energy fraction O.^, can be constrained 
by fitting the Hubble diagrams corresponding to such simulated distributions. It is the 
aim of this paper to simulate different GRB distributions and investigate their ability to 
determine the cosmological parameters Q. m and O.^. Both universes with and without a 
cosmological constant A will be considered. 



TO START: A KS TEST 



First, in order to show what we are aiming at, we performed a Kolmogorov-Smirnov 
(KS) test on two data sets made of 300 GRBs simulated in two different cosmological 
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FIGURE 1. Comparison of the cumulative fractions obtained with the parameter X = \ogd]^(Gpc) 
calculated for two data sets of 300 GRBs simulated in an universe with Q.,,, = 1 and £2a = (full line) and 
in one with fi„, = 0.3 and O.^ = 0.7 (dashed line). The corresponding probability that the two data sets 
are drawn from the same distribution is Q^s = 2.48 • 10~ 14 . 



models, one with Q. m = 1 and = and the other with Cl m = 0.3 and = 0.7, but 
both with a Hubble constant H$ = 65 km s Mpc 1 (as it will be assumed throughout 
the paper). We assume that GRBs are indeed standard candles with true prompt y- 
ray energy released, Ey, following a Gaussian distribution in its logarithm with mean 
/l = 50.7 (if Ey is expressed in erg units) and a = 0.3 (corresponding to a multiplicative 
factor of 2) ||2J], and that they are distributed in the universe according to the model 
of star formation rate Rsfi (z) reported in |0], which matches the log/V — logi 5 relation 
(GRB number counts vs. peak photon flux) obtained with BATSE data. Applying the KS 
test on the redshift distributions, we found that the probability that the two data sets are 
drawn from the same distributions is Qxs = 0.031, a "no man's land" value for this test. 
On the other hand, the application of the KS test on the parameter logd^(z), where di(z) 
is the luminosity distance, resulted in a significant probability Q^s ~ 10 -14 , which tells 
us that it is possible to discriminate between the two different cosmological models if a 
set of 300 GRB luminosity distances is known (see Figured)- 



DATA SET SIMULATIONS IN A A = COSMOLOGY 

We consider now a A = cosmology, in which the only contribution to the density 
parameter is given by Q. m . We assume for GRBs the same energy distribution as for the 
KS test. However, the assumed mean value is not relevant for our investigation, since it 
is the dispersion value that constrains the cosmological density parameter. 

The standard candle energy is related to the fluence of the burst f y = E r (l + 
z) I (4nd]^(z)) via the luminosity distance cIl(z). In order to have a linear propagation of 
errors throughout our simulations, we choose to construct with GRBs a Hubble diagram 
logd^—z, since the distribution of the parameter logJ^ is the same of that of log Ey, and 
therefore it is Gaussian. 

In order to study the ability of GRBs in probing the cosmological parameters as a 
function of their number, we have simulated different samples with Ngrb = 10, 30, 




TABLE 1. Mean values of the fitted cosmological density parameters £2,„ and £2a, of their error 
A£2 and their dispersion Sq obtained by fitting 10 2 GRB sample realizations with Ngrb distributed 
according to function Rsfi{z) of |3J] in an Einstein-de Sitter universe (ii„, = 1, left) and in a flat 
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10 


0.9983 


0.2997 


0.3097 


10 


0.3195 


0.6805 


0.1004 


0.1307 


30 


1.0158 


0.1895 


0.1993 


30 


0.2973 


0.7027 


0.0763 


0.0700 


100 


0.9937 


0.0993 


0.1108 


100 


0.3002 


0.6998 


0.0363 


0.0351 


300 


0.9959 


0.0599 


0.0629 


300 


0.3023 


0.6977 


0.0219 


0.0222 


1000 


1.0009 


0.0332 


0.0351 


1000 


0.3001 


0.6999 


0.0120 


0.0125 



100, 300 and 1000. Moreover, in order to be free from statistical fluctuations, we have 
performed 10 2 realizations of each of these samples. 

The simulation of a GRB consists of the random sampling of both the redshift z and 
the true y-ray energy released Ey, according to the respective adopted distributions. 
Given a cosmological model, from these coupled values we obtain the corresponding 
value for the parameter log J 2 , which we plot on the Hubble diagram as a function of 
Z. At this point we perform a x 2 minimization of the simulated data to see with which 
accuracy the fit reproduces the input cosmology. The measurement error on logJ^ is 
assumed to be o = 0.3. The mean results of our repeated fits in an Einstein-de Sitter 
universe (Q. m = 1) are reported in the left side of Table [T] 



DATA SET SIMULATIONS IN A A-DOMINATED COSMOLOGY 

We move now to a A-dominated cosmology, in which the contributions to the density 
parameter are given by the mass density, Q, m , and by the cosmological constant energy 
density, ^a. In light of the recent observations of the cosmic microwave background 
anisotropy [Q]], we put the prior of a flat universe Q. m + Q.\ = 1 . 

Again, in order to study the ability of GRBs in probing the cosmological parameters in 
a A-dominated universe, we have simulated 10 2 realizations of GRB samples with Ngrb 
= 10, 30, 100, 300 and 1000. The x 2 minimization of the resulting Hubble diagrams has 
been performed considering log d\ depending only on the fit parameter Q. m , i.e., using 
the relation ^a = 1 — Q m - The right side of Tabled reports the general results of our 
repeated fits for a flat cosmology with input values Q. m = 0.3 and ^a = 0.7 (which are 
those adopted in 

Focussing on the samples with Ngrb — 300, which represent the future data set 
expected from the Swift satellite experiment, Figure|2]shows one of the Hubble diagrams 
log d\ — z obtained with the simulations (left), together with the distribution of the best 
fit values of the matter density fraction Q. m for 10 3 sample realizations (right). 

Finally, we must remark that the analysis in ||2|] assumes of course a particular set 
of cosmological parameters to derive the standard y-ray energy of GRBs. To avoid a 
circular logic we should assume a candle calibration with a local sample of sources, a 
prospect which can now be considered possible in light of the discovery of the near GRB 
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FIGURE 2. Left: Hubble diagram logoff — z with data simulated for a sample of 300 GRBs in a fiat 
universe with density parameters Cl m = 0.3 and Q.a = 0.7. The solid curve shows the function logdf (z) in 
the assumed cosmology, while the dashed curves give the dispersion about the best fit parameter (upper 
curve corresponds to lower £l m ). Right: Histogram with the distribution of the best fit values of the matter 
density £l m for 10 3 realizations of a sample of 300 GRBs in a flat universe with density parameters 
Q. m = 0.3 and = 0.7. The distribution has a mean < Q. m >= 0.3001, a median £l m (med) = 0.3002, 
a dispersion Sq ih = 0.0228, and a kurtosis k^ m = 3.0993, to be compared with the value of a Gaussian 
distribution, i.e., 3. 



030329, with redshift as low as z = 0.1685. 



CONCLUSIONS 

We have simulated different samples of GRBs adopting y-ray energy and redshift distri- 
butions consistent with recent observational results, in order to investigate their ability 
to probe cosmological parameters such as the density fractions Q, m and Q.\. Our result is 
that in a A-dominated flat universe the accuracy in the determination of the matter den- 
sity Cl m is ~40% for a sample with Ngrb = 10 and an excellent ~4% for Ngrb = 1000. 
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